home *** CD-ROM | disk | FTP | other *** search
/ Amiga Format CD 46 / Amiga Format CD46 (1999-10-20)(Future Publishing)(GB)[!][issue 1999-12].iso / -in_the_mag- / reader_requests / scilab / tests / rtitr.tst < prev    next >
Text File  |  1999-09-16  |  2KB  |  68 lines
  1. //test de la fonction rtitr
  2. eps=100*%eps;
  3. //
  4. //siso
  5. //----
  6. //
  7. //causal
  8. //n1 scalaire
  9. n1=1;d1=poly([1 1],'s','c');       // yj=y(j-1)+u(j-1)
  10. r1=[0 1 0 1 0 1 0 1 0 1 0];
  11. r=rtitr(n1,d1,ones(1,10));if norm(r1-r)>eps then pause,end
  12. //redemarrage a chaud
  13. r=rtitr(n1,d1,ones(1,9),1,0);if norm(r1(2:11)-r)>eps then pause,end
  14. //n1 polynomial
  15. n1=poly(1,'s','c');
  16. r=rtitr(n1,d1,ones(1,10));if norm(r1-r)>eps then pause,end
  17. //
  18. r=rtitr(n1,d1,ones(1,9),1,0);if norm(r1(2:11)-r)>eps then pause,end
  19. //
  20. //non causal
  21. n2=poly([1 1 1],'s','c');d2=d1;    // yj=-y(j-1)+u(j-1)+u(j)+u(j+1)
  22. r2=[2 1 2 1 2 1 2 1 2];
  23. r=rtitr(n2,d2,ones(1,10));if norm(r-r2)>eps then pause,end
  24. //redemarrage a chaud
  25. r=rtitr(n2,d2,ones(1,9),1,2);if norm(r2(2:9)-r)>eps then pause,end
  26. //
  27. //mimo
  28. //----
  29. //
  30. //causal
  31. d1=d1*diag([1 0.5]);n1=[1 3 1;2 4 1];r1=[5;14]*r1;
  32. r=rtitr(n1,d1,ones(3,10));if norm(r1-r)>eps then pause,end
  33. //
  34. r=rtitr(n1,d1,ones(3,9),[1;1;1],[0;0]);if norm(r1(:,2:11)-r)>eps then pause,end
  35. //n1 polynomial
  36. n1(1,1)=poly(1,'s','c');
  37. r=rtitr(n1,d1,ones(3,10));if norm(r1-r)>eps then pause,end
  38. //
  39. r=rtitr(n1,d1,ones(3,9),[1;1;1],[0;0]);if norm(r1(:,2:11)-r)>eps then pause,end
  40. //non causal
  41. d2=d1;n2=n2*n1;r2=[5;14]*r2;
  42. r=rtitr(n2,d2,ones(3,10));if norm(r2-r)>eps then pause,end
  43. //
  44. r=rtitr(n2,d2,ones(3,9),[1;1;1],[10;28]);if norm(r2(:,2:9)-r)>eps then pause,end
  45.  
  46. //
  47.   a = [0.21 , 0.63 , 0.56 , 0.23 , 0.31
  48.        0.76 , 0.85 , 0.66 , 0.23 , 0.93
  49.        0 , 0.69 , 0.73 , 0.22 , 0.21
  50.        0.33 , 0.88 , 0.2 , 0.88 , 0.31
  51.        0.67 , 0.07 , 0.54 , 0.65 , 0.36];
  52.   b = [0.29 , 0.5 , 0.92
  53.        0.57 , 0.44 , 0.04
  54.        0.48 , 0.27 , 0.48
  55.        0.33 , 0.63 , 0.26
  56.        0.59 , 0.41 , 0.41];
  57.   c = [0.28 , 0.78 , 0.11 , 0.15 , 0.84
  58.        0.13 , 0.21 , 0.69 , 0.7 , 0.41];
  59.   d = [0.41 , 0.11 , 0.56
  60.        0.88 , 0.2 , 0.59];
  61. s=syslin('d',a,b,c,d);
  62. h=ss2tf(s);num=h(2);den=h(3);den=den(1,1)*eye(2);
  63. u=1;u(3,10)=0;
  64. r3=flts(u,h);
  65. //r=rtitr(num,den,u);if norm(r3-r)>1000*eps then pause,end
  66.  
  67.  
  68.